Carte inondation

0.1 Fond de carte

library(leaflet)
library(maptools)
library(rgeos)
library('rgdal')

alsace <- readShapePoly("data/inondation_shp/68_ATLAS_INONDATION_ZONE_INONDABLE",proj4string = CRS("+init=epsg:2154"))

alsace <- spTransform(alsace, "+init=epsg:4326")

leaflet() %>%addTiles()%>%
  addPolygons(data=alsace)

set.seed(1)
df <- sp::SpatialPointsDataFrame(
  cbind(
    7.3895+rnorm(10,0,0.1),  # lng
    48.0661+rnorm(10,0,0.1)  # lat
  ),
  data.frame(type = rep("a",10))
)


# sp::coordinates(df) <- ~Long+Lat
sp::proj4string(df) <- sp::proj4string(alsace)

points=sp::over(df,alsace)

p=df[is.na(points$IDENT)==FALSE,]
p

##             coordinates type
## 2  (7.407864, 48.10508)    a
## 3  (7.305937, 48.00398)    a
## 5  (7.422451, 48.17859)    a
## 7  (7.438243, 48.06448)    a
## 8  (7.463332, 48.16048)    a
## 9  (7.447078, 48.14822)    a
## 10 (7.358961, 48.12549)    a

leaflet() %>%addTiles()%>%
  addPolygons(data=alsace)%>%
  addMarkers(data=p)%>%
  addCircleMarkers(data=df,color="red")

dep_int <- readShapePoly("data/DEPARTEMENT/DEPARTEMENT",proj4string = CRS("+init=epsg:2154")) # ici nous avons du recharger le fond de carte departement avec le code 2154 qui signifie le format RGF93 / Lambert-93 de notre fond de carte de base, RGF93 :Réseau géodésique français légal pour la métropole depuis le décret du 26.12.2000
dep_int <- spTransform(dep_int, "+init=epsg:4326")

leaflet() %>%addTiles()%>%
  addPolygons(data=dep_int)

sp::over(df,dep_int)

##                   ID_GEOFLA CODE_DEPT  NOM_DEPT CODE_CHF NOM_CHF
## 1  DEPARTEM0000000000000005        68 HAUT-RHIN      066  COLMAR
## 2  DEPARTEM0000000000000005        68 HAUT-RHIN      066  COLMAR
## 3  DEPARTEM0000000000000005        68 HAUT-RHIN      066  COLMAR
## 4  DEPARTEM0000000000000005        68 HAUT-RHIN      066  COLMAR
## 5  DEPARTEM0000000000000005        68 HAUT-RHIN      066  COLMAR
## 6  DEPARTEM0000000000000005        68 HAUT-RHIN      066  COLMAR
## 7  DEPARTEM0000000000000005        68 HAUT-RHIN      066  COLMAR
## 8  DEPARTEM0000000000000005        68 HAUT-RHIN      066  COLMAR
## 9  DEPARTEM0000000000000005        68 HAUT-RHIN      066  COLMAR
## 10 DEPARTEM0000000000000005        68 HAUT-RHIN      066  COLMAR
##    X_CHF_LIEU Y_CHF_LIEU X_CENTROID Y_CENTROID CODE_REG
## 1     1024125    6784581    1016174    6763894       44
## 2     1024125    6784581    1016174    6763894       44
## 3     1024125    6784581    1016174    6763894       44
## 4     1024125    6784581    1016174    6763894       44
## 5     1024125    6784581    1016174    6763894       44
## 6     1024125    6784581    1016174    6763894       44
## 7     1024125    6784581    1016174    6763894       44
## 8     1024125    6784581    1016174    6763894       44
## 9     1024125    6784581    1016174    6763894       44
## 10    1024125    6784581    1016174    6763894       44
##                              NOM_REG
## 1  ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 2  ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 3  ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 4  ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 5  ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 6  ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 7  ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 8  ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 9  ALSACE-CHAMPAGNE-ARDENNE-LORRAINE
## 10 ALSACE-CHAMPAGNE-ARDENNE-LORRAINE